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The variational iteration method, a well-known method for solving functional equations, is em- 
ployed to solve a class of eighth-order boundary-value problems, which govern scientific and engi- 
neering experimentations. Some special cases of the mentioned equations are solved as examples to 
illustrate the ability and reliability of the method. The results reveal that the method is very effective 


and convenient. 
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1. Introduction 


The solution of nonlinear problems by analytical 
techniques is often rather difficult. Higher-order initial- 
boundary-value problems arise in many engineering 
applications [1-5]. Various numerical and analyti- 
cal methods were proposed to solve such problems 
(see [6 — 10] and references cited therein). Eighth-order 
differential equations govern physics of some hydrody- 
namic stability problems. When an infinite horizontal 
layer of fluid is heated from below and is subjected to 
rotation, instability is observed. When this instability 
starts as an overstability, it is modelled by an eighth- 
order ordinary differential equation [11]. Siddiqi and 
Twizell [12] used polynomial splines of degree six 
to solve some special cases of linear eighth-order 
boundary-value problems. However, the obtained re- 
sults were divergent at points adjacent to the boundary. 
Twizell and his coworkers [13,14] have also solved 
some other higher-order problems and encountered the 
same problems. Liu and Wu [15], utilizing the gener- 
alized differential quadrature rule, have solved similar 
problems. 

In this paper, a class of eighth-order boundary- 
value problems is considered. Four examples, solved 
in [12, 15], are dealt with again to obtain accurate re- 
sults in the entire domain via the variational iteration 
method. We consider the following class of eighth- 
order boundary-value problems: 


u®)(x)+0(x)u(x) =Via), a<x<b, 


with the boundary conditions 


u“)(a) = Ag, 
u)(b) = Bo, 


where u = u(x) and @(x) and y(x) are continuous func- 
tions defined on [a,b] and the constants A; and B; are 
finite real numbers. Furthermore, we will apply the 
variational iteration method on the following exam- 
ple [16]: 


O<x<l, (2) 


The article is organized as follows: Section 2 briefly 
reviews the mathematical basis of the variational it- 
eration method used for this study. Five illustrative 
examples are documented in Section 3. These exam- 
ples intuitively describe the ability and reliability of 
the method. A conclusion and future directions for re- 
search are all summarized in the last section. 


2. The Variational Iteration Method 


The basic concept of the variational iteration method 
is illustrated in [17-22], and applications of the 
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method can be found in [6—8, 23-31]. We consider 
the general nonlinear system 


Llu(x)] + Rlu(x)| + N[u()] = wa), 


where L = = R is a linear operator, N the nonlinear 
term and w(x) a given continuous function. According 
to the variational iteration method, we can construct a 
correction functional in the form 


x 
nA eos [ A.(s) [Latn(s) +Rai,(s) 

+ Niin(s) — W(s)]ds, 
where uo(x) is an initial approximation with possible 
unknowns, A is a Lagrange multiplier which can be 
identified optimally via the variational theory, the sub- 
script n denotes the nth approximation, and a, is con- 
sidered as a restricted variation, i.e., d#, = 0. 

For eighth-order boundary-value problems men- 
tioned above, according to the variational iteration 
method, the nonlinear terms have to be considered as 


restricted variation. So we derive a correction func- 
tional as follows: 


unsi(t) =tn(x) + f° Als) ®(s) 
+ 0(s)ii(s) — y(s)]ds. 


The stationary condition of the above correction func- 
tional can be expressed as 


1®)(s) =0, 
1- A (5) sax = 0, 
AO) [eae 0; FH 152500446, 
A(s)|s=x = 0. 
The Lagrange multiplier, therefore, can be identified as 
A= as ae 


As aresult we obtain the iteration formula 
1 x 
tnsi(a) = tls) +27 [(s—x)" WC) 
+ 6(s)u(s) — y(s)]ds. 


In the above-mentioned iteration formula, we begin 
with an arbitrary initial approximation: 


(3) 


uo(x) =a+bx+cx? +dx3 + kxt + fx + gx°+ hx’, 


where a, b, c, d,k, f, g and h are constants to be deter- 
mined. 
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3. Applications 


In this section, we present five examples to show the 
efficiency and high accuracy of the present method. 


3.1. Example 1 


Consider (1) with [a,b] = [0,1], (x) = x, 
w(x) = —(48 + 15x+.x3)e*, and the boundary condi- 
tions [12, 15] 

Ao = 0, A2 = 0, 
Bo = 0, Bo => —4e, 


As = —8, 
Ba = —16e, 


Ag = —24, 
(4) 
Be = —36e. 


According to (3) we have the following iteration for- 
mula: 


Unt1(X) = Un(x) + = [o—ay"tu (s) 


+ sun(s) + (48 + 15s+53)e*]ds. 


(5) 


Now, we begin with the arbitrary initial approximation 
uo(x) = a+bx+ cx? + dx? + kxt + fxr + gx®+ hx’, 


where a, b, c, d, k, f, g and h are constants to be deter- 
mined. 
By the variational iteration formula (5), we have 


uy (x) = ug(x) + - [o —x)"[u)(s) + suo(s) 
+(48 + 15s+s3)e']ds 


79 
= a—792+ (b—561)x+ (c — 189)x? + («-4)x 
1D 19\ , 1: 
oy («- 5 )a +(r-3)s m G oa) 
3 7 ax? bx cx! 
Riot yee eens fore So eae ey 
560 362880 1814400 6652800 
dx'2 kx}3 al gx 
~ 19958400 51891840 121080960 259459200 
hx'6 
S006 = 3 se ae ae, 
sisciging ee ee 


Incorporating the boundary conditions (4) into u(x) 
yields a linear system with 8 equations and 8 variables. 
Solving this linear system simultaneously, we have a = 
0,b=1c=0d=53, k=, f=y.¢ = 5 and 
h= ah. Thus we obtain the first-order approximate 
solution 


Re. ape Sy x 
2 3 8 
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Table 1. Comparison of the first-order approximate solution, 
ui (x), with the exact one. 


x Exact First-order Absolute 
solution approximate solution error 

0 0 0 0 

0.1 0.9946538262 0.9946538265 <10-° 
0.2 0.1954244413 0.1954244445 <10-8 
0.3 0.2834703497 0.2834704312 < 10-6 
0.4 0.3580379275 0.3580387556 < 10-5 
0.5 0.4121803178 0.4121853299 < 10-5 
0.6 0.4373085120 0.4373304000 < 10-4 
0.7 0.4228880685 0.4229643785 < 10-3 
0.8 0.3560865485 0.3563121778 < 10-3 
0.9 0.2213642800 0.2219525438 < 10-3 
1 0 0.0013888888 <10-? 


0 0.2 0.4 0.6 0.8 i 
x 
Fig. 1. Comparison of the first-order approximate solution, 
u(x), with the exact one. Solid line, exact solution; symbols, 
approximate solution. 


As the same, we find u(x) as follows: 


up (x) =m) ['(o—2)tul? (6) +5069) 
+(48 + 155+ s3)e']ds 


- a x4 x x x! x8 x? 
~*~3~ 37 8 30 144 840 5760 
x10 xl x2 x)3 
~ 45360 403200 3991680 43545600 
xl4 x15 


ee eS: 16 
518918400 6706022400 ae 


This gives the solution in a closed form by x(1 — x)e*. 
The first-order approximate solution is of remarkable 
accuracy. Comparison of the first-order approximate 
solution, u(x), with the exact one is shown in Table | 
and Figure 1. 


TA7 
3.2. Example 2 
Consider (1) with [a,b] = [—1,1], (x) = —x and 


w(x) = —(55 + 17x+.x7 —x3)e*, with the boundary 
conditions [12, 15] 


Aimy Aa 2/0, AeA Re Ae 181, a, 
Bo =0, By =—6e, Ba=—20e, Bo = —42e. 


According to (3) we have the following iteration for- 
mula: 


1 x 
tnsi(s)=tlt)-t ae f° (s— a)" ult) (5) — sul) 
+ (55+17s+5s?—s?)e'|ds. 


(7) 


Now, we begin with the arbitrary initial approximation 
uo(x) =a+bx+ cx? + dx? + kx? + fxr + gx°+ hx’, 


where a, b, c, d, k, f, g and h are constants to be de- 
termined. By the variational iteration formula (7), we 
have 


uy (x) = w(x) + af — x) [uf (s) — suo(s) 
+(55+ 17s +s" —s?)e5]ds 


231 
=711+a+ (496+ b)x+ (c+ >) x 4+ (d+35)x23 


127), 4 19\-< vie ee 
+ (or 2) te (42) 4 (e428) 


h 23 ae ax 4 bx10 & cx 
if pte | ge i ct 
2520 362880 1814400 6652800 
dx'2 kx}3 |e ge 
+ ——— 4+ ——_—— 4 oeoioe— Ff 
19958400 51891840 121080960 259459200 
hx'6 
518918400 


— 25e%x* +215e%x —71le*+e*x?. 


Incorporating the boundary conditions (6) into 
ui(x) yields a linear system with 8 equations and 
8 variables. Solving this linear system simultaneously, 
Welwed=1,b=1¢=—/d=2,h= 4,5 
f= = g= = and h = sa: Thus we obtain the 


first-order approximate solution 


ee Se Tl 4 19 2 


=] _ — 
2 a al a a Ya 
29, 41 , 
~ 720" ~ 5040" 
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Table 2. Comparison of the first-order approximate solution, 


ui (x), with the exact one. 


x Exact First-order Absolute 
solution approximate solution error 
-1 0 0.1190476190 < 107? 
—0.8 0.1617584271 0.1619634388 <10-3 
—0.6 0.3512394471 0.3512605257 <10-+ 
—0.4 0.5630688386 0.5630696838 < 10-5 
—0.2 0.78598 15230 0.78598 15264 <10-8 
0 1 1 0 
0.2 1.172546648 1.172546651 <1077 
0.4 1.253132746 1.253133695 < 10-5 
0.6 1.166156032 1.166181074 <10-+ 
0.8 0.8011947341 0.8014527387 < 10-3 
1 0) 0.1587301587 <10~? 


Fig. 2. Comparison of the first-order approximate solution, 
u(x), with the exact one. Solid line, exact solution; symbols, 
approximate solution. 


As the same, we find u2(x) as follows: 


up (x) =m) [(o—x) ul? (6) —s0() 
+(55+17s+5?—s3)e']ds 


2 6 24 120 720 5040 
lig. Fig? 89x10 109x!! 
~ 8064 362880 3628800 39916800 
1314 31x" 18ia" 
~ 479001600 1245404160 87178291200 
194° ie 
~ Trse7oass000 OW 


This gives the solution in a closed form by (1 —x7)e*. 
The first-order approximate solution is of remarkable 
accuracy. Comparison of the first-order approximate 
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solution, u;(x), with the exact one is shown in Table 2 
and Figure 2. 


3.3. Example 3 


Consider (1) with [a,b] = [—1,1], (x) = —1 and 
w(x) = —8[2x cos(x) + 7sin(x)], with the boundary 
conditions [12, 15] 

Ag =0, Az = —4cos(1) —2sin(1), 
Ag = 8cos(1) + 12sin(1), 
Ag = —12cos(1) — 30sin(1), 


Bo=0, Bz =4cos(1)+2sin(1), 8) 
Bg = —8cos(1) — 12sin(1), 
Bo = 12cos(1) + 30sin(1). 
According to (3) we have the iteration formula 
Lf T,,,(8) 
tnt) =tn(a)+a, f(s—a){ul(s) mals) og 


+ 8[2scos(s)-+7 sin(s)]}ds. 
Now, we begin with the arbitrary initial approximation 
uo(x) = a+bx+ cx? + dx? + kx* + fx + gx®+ hx’, 


where a, b, c, d, k, f, g and h are constants to be de- 
termined. By the variational iteration formula (9), we 
have 


= [= 2) {ul (6) — wos) 
+8[2scos(s) + 7sin(s)] }ds 


1 
=a+t (b—56)x+ ex? + (d+4)xP+ kext + (r+ a) a 
1 ax® 
6 a 
4 | pores 
grt ( a) * ¥ 79320 * 362880 
cx!0 dx'! kx!2 fx? 
+ 4} 
1814400 | 6652800 | 19958400 ' 51891840 
gx'4 hx'5 
———————— + a, 
121080960 | 259459200 


bx? 


— 16cos(x)x + 72 sin(x). 


Incorporating the boundary conditions (8) into u(x) 
yields a linear system with 8 equations and 8 variables. 
Solving this linear system simultaneously, we have a = 
0,b=-1,c=0,d=%,k=0, f= 7, g =0 and 
h= a Thus we obtain the first-order approximate 
solution 


7x 7 


6 40° 


5 
+ 50407 
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Table 3. Comparison of the first-order approximate solution, 
ui (x), with the exact one. 


x Exact First-order Absolute 
solution approximate solution error 
-1 0 —0.1984126984 < 10-3 
—0.8 0.2582481927 0.2582214298 < 107 
—0.6 0.3613711830 0.3613691657 <10-+ 
—0.4 0.3271114075 0.3271113549 <10~6 
—0.2 0.1907225576 0.1907225575 <10~°° 
0 0 0 0 
0.2 —0.1907225576 —0.1907225575 <10~-° 
0.4 —0.3271114075 —0.3271113549 < 10-6 
0.6 —0.3613711830 —0.3613691657 <10-+ 
0.8 —0.2582481927 —0.2582214298 <10-+ 
1 0 0.1984126984 < 1073 


Fig. 3. Comparison of the first-order approximate solution, 
ui (x), with the exact one. Solid line, exact solution; symbols, 
approximate solution. 


Similarly we can find u(x) as follows: 


u2(x) = u(x) + = [ou —uy(s) 
+8[2scos(s) +7 sin(s)|}ds 


ie T° 43x’ 73x? 37x!! 
=> —X —_ - —_ — —_—_—_ 
6 40 | 5040 362880 13305600 
isqa* ia 
SG 
@227020800 | O""”? 


This gives the solution in a closed form by (x? — 1) 
sin(x). The first-order approximate solution is of re- 
markable accuracy. Comparison of the first-order ap- 
proximate solution, u(x), with the exact one is shown 
in Table 3 and Figure 3. 


3.4. Example 4 


Consider (1) with [a,b] = [—1,1], (x) = —1 and 
w(x) = 8[2x sin(x) —7cos(x)], with the boundary con- 
ditions [12, 15] 

Ap =0, Az =—4sin(1)+2cos(1), 
Ag = 8sin(1) — 12cos(1), 

Ag = —12sin(1) + 30cos(1), 
Bo=0, Bo = —4sin(1)+2cos(1), 
By = 8sin(1) — 12cos(1), 

Bo = —12sin(1) + 30cos(1). 


(10) 


According to (3) we have the iteration formula 


unei()=unlat ze [(s—2)" (ul) —Un(s) (11) 


— 8[2ssin(s)—7cos(s)|}ds. 
Now, we begin with the arbitrary initial approximation 
ug(x) = a+bx+ cx? + dx? + kx! + fo +.9x°+ hex’, 


where a, b, c, d, k, f, g and h are constants to be de- 
termined. By the variational iteration formula (11), we 
have 


uy(x) = u(x) + = [0 tubo) —uo(s 
—8[2ssin(s) — 7cos(s) 


) 
I}as 
1 
=—72+a+bx+(c+20)x?+dx + («- :) - 


ax® bx? 


1 
5 6 7 
ae — h 
fx ( =) af “ * 40320 ' 362880 


ky!2 Fe 
4+} 
19958400 51891840 


cx!0 dx!! 


*Tg14400 * 6652800 


exl4 nels 


* 731080960 | 259459200 


Incorporating the boundary conditions (10) into u; (x) 
yields a linear system with 8 equations and 8 variables. 
Solving this linear system simultaneously, we have a = 
=Lb=0,c=2,d=0k= 2,7 =0,2=4, adh= 
0. Thus we obtain the first-order approximate solution 


+ 16sin(x)x +72 cos(x). 


2-4 1 
uy(x) = — Pee an ces 


2 24 720 


Similarly we can find u2(x) as follows: 


uz(x) = uy(x) + = [ou —uy(s) 
—8[2ssin(s) — 7cos(s)]}ds 
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Table 4. Comparison of the first-order approximate solution, 
ui (x), with the exact one. 


x Exact First-order Absolute 
solution approximate solution error 
-1 0 —0.1388888889 < 107? 
—0.8 —0.2508144153 —0.2505799111 < 10-3 
—0.6 —0.5282147935 —0.5281912000 < 10-3 
—0.4 —0.7736912350 —0.7736903111 <10-5 
—0.2 —0.9408639147 —0.9408639111 <10-8 
0 —1 -1 0 
0.2 —0.9408639147 —0.9408639111 <10-8 
0.4 —0.7736912350 —0.7736903111 <10>5 
0.6 —0.5282147935 —0.5281912000 < 10-3 
0.8 —0.2508144153 —0.2505799111 < 10-3 
1 0 —0.1388888889 < 107? 


x 


Fig. 4. Comparison of the first-order approximate solution, 
u1(x), with the exact one. Solid line, exact solution; symbols, 
approximate solution. 


= 143% 13x4 31x® 1938 a Tax 
7 2 24 720 13440 ° 518400 
1 12 1 14 
ee a ne i Ole. O(«!®). 
68428800 — 29059430400 


This gives the solution in a closed form by (x? — 1) 
cos(x). The first-order approximate solution is of re- 
markable accuracy. Comparison of the first-order ap- 
proximate solution, u1(x), with the exact one is shown 
in Table 4 and Figure 4. 


3.5. Example 5 


In order to illustrate the effectiveness of the vari- 
ational iteration method, we consider the nonlinear 
boundary-value problem [16] 

u® (x) =e *u"(x), 


O0<x<1l, (12) 
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Table 5. Comparison of the first-order approximate solution, 
u(x), with the exact one. 


x Exact First-order Absolute 
solution approximate solution error 
0.1 1.105170918 1.105170918 <10-° 
0.2 1.221402758 1.221402759 <10-8 
0.3 1.349858808 1.349858805 <10-8 
0.4 1.491824698 1.491824681 <10-7 
0.5 1.648721271 1.648721168 < 10-6 
0.6 1.822118800 1.822118354 < 10-6 
0.7 2.013752707 2.013751158 <10> 
0.8 2.225540928 2.225536366 <10-4 
0.9 2.459603111 2.459591262 <10-3 
1 e 2.718253968 < 1074 


ie) 0.2 0.4 0.6 0.8 1 


Fig. 5. Comparison of the first-order approximate solution, 
u(x), with the exact one. Solid line, exact solution; symbols, 
approximate solution. 


with the boundary conditions 

0)=1, 
2 4 6 (13) 

u(1) =e, u?)(1) =e, u(1) =e, uw (1) =e. 


According to the variational iteration method we have 
the following iteration formula: 


tne1() uals) +a [ (s=2)"tul® (9) eu] 
(14) 


Now, we begin with the arbitrary initial approximation 
uo(x) =a+bx+cx?+ dx? + kx* + fx + gx®+ hx’, 


where a, b, c, d, k, f, g and h are constants to be de- 
termined. By the variational iteration formula (14), we 
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obtain u(x). Because of exhaustive computations, we 
do not write uj (x). Incorporating the boundary condi- 
tions (13) into u(x) yields a linear system with 8 equa- 
tions and 8 variables. Solving this linear system simul- 
taneously, we havea=1,b=1,c 5. d é k a 
f= ae g= aT andh= un Thus we obtain the first- 
order approximate solution 


Cee ee ee le 
uy(x) = x+—-+—+-—4+—x 
: g." § D4. 190 
1 6 1 7 
+ 590% + 50407 ° 


As the same, we find w(x) as follows: 


+ = |G =x) u® (s) —e~Sui(s)]ds 


5 


u2(x) = uy (x) 


ae al 1 1 1 
= Bp Bg. toa No Sy te 
TR gh TEN gat Gogh agg 
1 4 1 of 9 
a ae BO (98) 
+ sag" + qo300" TO) 
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4. Conclusions 


We have studied some eighth-order boundary-value 
problems using the variational iteration method. The 
initial approximation was selected as a polynomial 
with unknown constants, which was determined by 
considering the boundary conditions. The results re- 
vealed that the method is remarkably effective. This 
method is a very promoting one, which will be cer- 
tainly found widely applied. 
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